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ABSTRACT 

A natural architecture for nanoscale quantum computation is that of a quantum cellular au- 
tomaton. Motivated by this observation, in this paper we begin an investigation of exactly 
unitary cellular automata. After proving that there can be no nontrivial, homogeneous, 
local, unitary, scalar cellular automaton in one dimension, we weaken the homogeneity con- 
dition and show that there are nontrivial, exactly unitary, partitioning cellular automata. 
We find a one parameter family of evolution rules which are best interpreted as those 
for a one particle quantum automaton. This model is naturally reformulated as a two 
component cellular automaton which we demonstrate to limit to the Dirac equation. We 
describe two generalizations of this automaton, the second of which, to multiple interacting 
particles, is the correct definition of a quantum lattice gas. 
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1. Introduction 

The realization that incentives to develop smaller and faster computers will eventually 
drive the devices from which they are constructed into the quantum regime motivated 
research into quantum mechanical limitations on deterministic computation as early as 
the 1970s [1]. The subsequent conceptions of universal quantum simulator by Feynman 
[2] and quantum Turing machine by Deutsch [3] initiated a series of investigations [4] into 
how aspects of quantum mechanics, specifically superposition and interference, might be 
exploited for computational purposes. Shor's remarkable discovery of a polynomial time 
quantum algorithm for factorization [5] (and the centrality of the factoring problem to 
modern cryptography [6]) has led to redoubled interest in the design and construction of 
quantum computational nanodevices [7]. It should be emphasized that the goal here is 
a computational device which will run quantum algorithms, not a quantum device which 
will run deterministic [8,9] or probabilistic [10] algorithms. 

For a variety of reasons — the wire and gain problems, and the pragmatic observation 
that an array of simple devices is often easier to design and build than a single, more com- 
plicated device — it seems likely that massive parallelism will optimize nanoscale computer 
architecture [11]. In this paradigm, a quantum computer is a quantum cellular automaton 
(QCA): the state of each simple device (cell) in the array depends on the states of the 
cells in some local neighborhood at the previous timestep. Unlike the original cellular 
automaton (CA) models for parallel computation of von Neumann and Ulam [12], where 
this dependence is deterministic or probabilistic, here the dependence is quantum mechan- 
ical: There is a (complex) probability amplitude for the transition to each possible state, 
subject to the condition that the evolution be unitary, so that the total probability — the 
sum of the norm squared of the amplitude of each configuration — is always one. 

QCA, therefore, provide a laboratory for analyzing both potential quantum computer 
architectures and algorithms; this is the motivation for initiating our study of them in this 
paper. Computation motivated study of QCA seems to have originated with the interest- 
ing work of Grossing and Zeilinger [13]. Their models, however, are only approximately 
quantum mechanical because, they argue, "except for the trivial case, strictly local, uni- 
tary evolution of the whole QCA array is impossible" [14] . Consequently, they study CA 
whose evolution is both nonunitary [13] and, although 'probability' preserving, nonlinear 
[14] . In Section 2 we begin by proving the following precise version of their claim: 

NO-GO Lemma. In one dimension there exists no nontrivial, homogeneous, local, scalar 
QCA. More explicitly, every band r -diagonal unitary matrix which commutes with the 
1-step translation matrix is also a translation matrix, times a phase. 

and then continue by showing that even a slight weakening of the homogeneity /translation 
invariance condition allows nontrivial unitary evolution. Thus we will reserve the adjective 

'quantum' for CA with exactly unitary, nontrivial, local evolution, in contrast to the usage 
by Grossing, Zeilinger, et al. [13,14] and Lent and Tougaw [9]. These are the QCA which 
should best model truly quantum parallel architectures. In this paper we consider only the 
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one dimensional situation; but this is not irrelevant for computational complexity issues 
since reversible CA capable of universal computation exist in one dimension [15], and 
reversible deterministic CA are, of course, unitary. 

Section 3 contains output from several simulations of one dimensional QCA. Simu- 
lation on a deterministic computer must inevitably be slow, but this is acceptable since 
our goal here is understanding rather than the solution of any specific problem. It is easy 
to see qualitative differences from the simulations of Grossing, Zeilinger, et al. [13,14]; in 
particular, our simulations display particle-like features. 

The latter observation motivates our reinterpretation in Section 4 of the evolution 
rule of such a QCA as the scattering rule for a quantum particle automaton. In this form 
the model is equivalent to Feynman's path integral formulation for a Dirac particle [16]; it 
is straightforward to solve the model exactly and to give a physical interpretation to the 
parameter in the scattering rule. 

A lattice gas [17] formulation of this quantum particle automaton would consist of an 
array of nodes occupied by left and/or right moving particles which jump to the next node 
at the next timestep. Since there is only one particle in the model the amplitudes for left 
and right moving particles at a given node may be combined into a two component field 
as in the one dimensional Dirac equation in the chiral representation. Formulated as such 
a two component QCA, our quantum particle automaton is unitarily equivalent to the 
quantum lattice Boltzmann equation of Succi and Benzi [18] and to the one dimensional 
version of Bialynicki-Birula's QCA [19]. Since a two component QCA evades the conclusion 
of the No-go Lemma, in Section 5 we find the most general homogeneous/translation 
invariant unitary evolution rules for a neighborhood of radius one. The lattice gas paradigm 
motivates another generalization, however, to multiple (interacting) quantum particles. 
This new system, described in Section 6, is a quantum lattice gas automaton (QLGA) and 
may be expected to be relevant for modelling parallel quantum architectures constructed 
from few electron devices. 

As Landauer [20] and others [21] have emphasized, achieving practical quantum com- 
putation will be difficult for a variety of reasons. The most fundamental problem is de- 
coherence [22,23], which destroys the delicate interference phenomena on which quantum 
algorithms such as Shor's [5] depend. In Section 7 we observe that this is among the im- 
portant issues in quantum computation which can be investigated using these models and 
also remark on some related research directions. 

2. Quantum cellular automata 

A CA consists of a lattice L of cells together with a Geld (p : N x L ^ S, where N denotes 
the non-negative integers labelling timesteps and S is the set of possible states in which the 
field is valued. The field evolves according to some local rule, i.e., (j) satisfies a recurrence 
relation of the form 

^t+i{x) = f{M^ + e) I e e E{t,x)), (1) 
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where E{t,x) is a set of lattice vectors defining local neighborhoods for the automaton. 

In the Schrodinger picture of quantum mechanics the state of a system at time t is 
a state vector in some Hilbert space (see, e.g., [24]). The state vector evolves locally and 
unitarily, i.e., 

4>i+i = Uct>t, (2) 



where U is a. unitary matrix (more precisely, operator: UU^ = I = U^U). Thus, if the 
configuration space is discrete, and the Hilbert space has a computational basis \x), x & L, 
it is natural to try to construct a QCA model where </>t(x) is the (complex scalar) coefficient 
of \x) in (pt and the local evolution rule (1) is the unitary evolution rule (2) in this basis. 
Notice that this forces the QCA to be additive in Wolfram's terminology [25], i.e., (1) 
becomes 

4>t+i{x) = w{t,x + e)(j)t{x + e), (3) 

eeE{t,x) 



where the coefficients w{t,x + e) are constrained by the unitarity condition. If both E{t, x) 
and w{t,x + e) are independent of t and x, the CA is homogeneous. 

In one dimension (pt can be written as a column vector with ordered entries 4>t{x), in 
which case locality corresponds to U being band diagonal. If, furthermore, the QCA is 
homogeneous, U must be invariant under the action of the translation operator T on the 
lattice — the permutation matrix with Is on the subdiagonal — and (2) becomes 



0t+i(-l) 
0t+i(O) 
0t+i(+l) 



V 



/ 



W-r • ■ ■ 

W-r . . . W+r 

W-r . . . W+r 



( 



\ 



</'t(-l) 

0t(O) 
0t(+l) 



V 



(4) 



where w{t,x -\- e) = We- (With periodic boundary conditions, of course, the top and 
bottom r rows will wrap around.) Grossing and Zeilinger consider the case r = 1 but find 
no nontrivial unitary matrix U . More generally, 

NO-GO Lemma. In one dimension there exists no nontrivial, homogeneous, local, scalar 
QCA. More explicitly, every band r-diagonal unitary matrix U which commutes with the 
1-step translation matrix T is also a translation matrix for some k E times a phase. 

Proof. When r — 0, U — wqI and unitarity implies |woP = 1 so C/ is the 0-step 
translation matrix (the identity) times a phase. Assume that the statement is true 
for r — 1. For a 1-step translation invariant band r-diagonal matrix, unitarity implies: 
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W-rW-r + ly-r+ll^-r+l + " H Wj.-\Wr-l + Wj-Wj. = 1 

W-r+lW-r + • 1- Wr-lWr-2 + WrWr-1 = 

Wr-lW-r + WrW'-r+l = (5r-l) 
WrW-r = 0, (5r) 

together with the complex conjugate equations. By equation (Sj.), at least one of Wr and 
W-r vanishes, say Wr- Then we may assume W-r 7^ 0; otherwise the conclusion follows 
immediately from the inductive hypothesis. But then equation (S^-i) forces Wr-i = 0; 
. . .; equation (5_r+i) forces W-r+i = 0; and equation (5-^) becomes I^W-rP = I7 ^-c, the 
only nonzero weight, W-ri is a phase. Then U — W-rT^ ■, where is the r-step translation 
matrix. ■ 



(5-r) 



CA which evolve simply by translation are not very interesting, so Grossing and 
Zeilinger relax the unitarity constraint, setting 

W-i = iS, Wo = 1, w+i = iS, 

so that the evolution is only approximately unitary, with errors of 0(|5p) accumulating 
at each timestep [13]. Since the evolution is nonunitary, quantum probability is not pre- 
served. Instead, all the amplitudes are normalized by an overall factor at each step to 
make = 1; this may be thought of as a nonlocal (and non-quantum mechani- 

cal) aspect of the evolution [14]. 

Our choice, instead, is to weaken the homogeneity condition but to insist still on 
exactly unitary, local evolution, thus maintaining consistency with standard quantum me- 
chanics. The evolution matrix in (4) is 1-step translation invariant, i.e., TUT~^ = U; a, 
natural way to relax this constraint is to require T'^UT~'^ = U, so that the evolution is 
only 2-step translation invariant. For r = 1 we have then: 



U 



a b c 
d e f 
a b c 
d e f 



\ 



and unitarity requires: 



ad + bb + cc 
bd + ce 
cd 



1 dd + ee + ff = 1 
ed + fb = 

fd = 0, 



(6) 
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together with the complex conjugate equations. 

There are two types of solutions to these equations. The uninteresting ones have only 
a and d (or c and /) nonzero, both with norm 1; in this case the evolution is by translation 
and multiplication by alternating phases. The interesting solutions to (6) have c = d = 
(or a = / = 0) and the matrix 




unitary; in this case U is block diagonal, each block acting only on a pair of adjacent cells. 

Evolving by U at each timestep would partition the CA into a set of noninteracting 
systems each comprised of a pair of adjacent cells. Instead, we evolve by U and by TUT~^ 
at alternating timesteps, changing the pairing and allowing propagation. Such alternat- 
ing evolution has been referred to as a staggered rule for a checkerboard model in the 
context of probabilistic CA and two dimensional statistical mechanics models [26], and as 
a partitioning CA in the context of reversible CA [27]. Since unitary evolution includes 
deterministic reversible evolution, the No-go Lemma applies in the latter context in one 
dimension (of course, the overall phase referred to in the statement of the lemma must be 
1), and it is natural that triviality has been evaded there in the same way. 



3. Simulations 

The matrix S E U{2) may be parameterized as 

Q _ f e*" sin 6 e*^ cos 6 \ 
\e'^cose e*'^sin6' J ' 

with a — /3 — 7 + 5 = n (mod 27r). Rather than continuing in full generality, wc shall impose 
parity invariancc, which forces b = e and a = f (or b = e and c = d) in (7). Dividing out 
an overall phase which is unobservable {i.e., has no effect on probabilities), the cell pair 
evolution matrix S becomes 

g _ f i sin 9 cos 9 \ 

\ cos 9 isva.9 J ' ^ ' 



Simulation of a single timestep of this QCA is achieved by a series of matrix multipli- 
cations: 

f (t)t+i{x ^ f isin9 cos9\ f (t)t{x-l)\ 
\ (t>t+i{x) J \cos9 ism9j\ ^t{x) 

for each x = t-^l (mod 2), say. Although these matrix multiplications commute since they 
evolve disjoint pairs of cells, sequential simulation of this intrinsically parallel evolution 
will be slow, as it always is for CA. Furthermore, since the field values and the evolution 
parameters are complex numbers, the usual efficiencies of bitwise computation are un- 
available. Nevertheless, simulation of small systems is easy and informative; an additional, 
more serious difficulty will not appear until the simulations of QLGA in Section 5. 
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have been imposed. 



Figure 1 shows simulations of our QCA for 9 = 7r/6 and 9 = 7r/3, starting from the 
same random initial condition. The darkness of each cell is (positively) proportional to its 
probability, i.e., the norm squared of its amplitude (field value), and it appears that there 
is a probability fiow in each case, slower for 9 = 7r/3. This is expected since for 9 = 7r/2, 
U is proportional to the identity (the overall factor of i is unobservable) and there is no 
flow, while for 9 = 0, S simply interchanges the states of adjacent cells so that probability 
propagates with speed 1 in lattice units. 

Starting with the simplest symmetric initial condition demonstrates the probability 
propagation more clearly. In the simulation shown in Figure 2, the initial condition is 

^(^) = |l/V2 if a; e {15, 16}; 
i otherwise. 

Here 9 = tt/A. Measuring the locations of the two peaks of the probability distribution at 
successive timesteps indicates that the propagation speed is approximately 2/3 in lattice 
units. 

This is consistent with the evolution of the even simpler initial condition shown in Fig- 
ure 3. Here the only nonzero initial value is at a; = 0. Again 9 = and the propagation 
speed still appears to be approximately 2/3 in lattice units. In both of these simulations 
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there are peaks in the probabihty distribution which remain substantially localized for 
the duration of the evolution shown. This behavior, particularly in the symmetric sim- 
ulation shown in Figure 2, should be contrasted with the results of Grossing, Zeilinger, 
et al. [13,14], which demonstrate quite different qualitative features: In their simulations 
macroscopic patterns develop and there is nothing which has the particle-like appearance 
of the persistent localized peaks in these probability distributions. 



4. Quantum particle automata 



The underlying particle nature of this al- 
ternating/partitioning QCA can be seen 
on the trajectory lattice dual to the space- 
time lattice of the CA. As shown in Fig- 
ure 4, each adjacent pair of cells acted on 
by S has a dual pair of spacetime edges 
which intersect and then continue along 
the same trajectories at the next timestep, 
becoming the spacetime edges dual to the 
same pair of cells. On the spacetime tra- 




Figure 4. The spacetime CA lattice is the rect- 
angular array of square cells; the dual particle 
trajectory lattice consists of the diagonal edges. 



jectory lattice the alternating action of what should now be called the scattering matrix S 
is automatic; the values of (f) are attached to the edges of the trajectory lattice and undergo 
a unitary transformation by S at each vertex. 

The value of on a left/right pointing edge of the trajectory lattice should be inter- 
preted as the amplitude for a left/right moving particle being there. That is, equation (9) 
means that if there is a right moving particle at time t, at time t -|- 1 it either becomes a 
left moving particle with amplitude isin^ or continues moving right with amplitude cos^. 
In particular, if there is some nonzero amplitude on each of two intersecting right and left 
going edges of the trajectory lattice, they may be evolved independently and then added. 
There is no interaction, so this is effectively a model of a single particle; the QCA is a 
quantum particle automaton, (ptix) is the amplitude for the particle being in state |x) 
during the time interval {t,t + 1), where, if x takes integer values on the vertical lines in 
Figure 4, 

{left moving from x + 1 iit^x (mod 2) 
right moving from x iit = x (mod 2). 

That is, the computational basis \x) is a set of eigenstates of an operator measuring more 
than just position. Projecting onto the position subspace adds (incoherently, i.e., the 
probabilities add) the amplitudes for left and right movers existing at (t, x), for t = x 
(mod 2), and gives a much clearer picture of the evolution. Figures 5 and 6 show this 
representation of the simulations of Figures 2 and 3, respectively. Since the evolution lies 
within the 'lightcone', i.e., the propagation speed is less that 1 in lattice units. Figures 5 
and 6 are in lightcone coordinates: 

:= {t+{x-xo))/2 
{t-{x-xo))/2, 



\x) 



(10) 



u 
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Figure 5. Position measurements in liglitcone 
coordinates for the QCA of Figure 2. The prob- 
ability 1 measurement at (0, 0) has been clipped. 



Figure 6. The same representation of the QCA 
of Figure 3. The first few (large) probabilities in 
the evolution have been clipped. 



where xq implements a translation of the origin: by 16 in Figure 5 and by in Figure 6. In 
both figures only the spacetime region covered by < w, f < 16 is shown, the origin is at 
the rear corner, and time runs forward to the front corner. The parity dependence (visible 
in the patterns of alternating dark and light cells) in Figures 2 and 3 has been smoothed out 
in these so that both the probable particle trajectories as well as their wave-like character 
are clearly visible. 

That the solution is as smooth as shown here suggests that the quantum particle 
automaton may be a discrete approximation to a continuum system, possibly with an 
exact solution. This is, in fact, the case. Although we have been led to it simply by the 
assumptions of discreteness, locality, and unitarity (and parity invariance), we will show 
that this quantum particle automaton is a discrete approximation to the Feynman path 
integral for a Dirac particle in one dimension [16]. Furthermore, it is exactly solvable, even 
without going to the continuum limit. 



Consider a right moving particle at 
(tjx) = (0,0), i.e., (/>o(0) = 1. The ampli- 
tude 4>t{x) for the particle being in state 
I a;) at time t is defined to be the propaga- 
tor K{t, x;0,0). This is exactly the (co- 
herent) sum of the amplitudes of all pos- 
sible paths in the trajectory lattice from 
(0, 0) to {t, x) — the sum-over-histories, or 
Feynman path sum. Figure 7 shows two 
typical paths a right moving particle from 
(0, 0) might take, arriving in left moving 
state \x — 1) or right moving state \x) 
just after time t {= x (mod 2)). Each 
path consists of u steps right and v steps 
left, followed by the terminal step, and is 




Figure 7. A typical initially right moving path 

from (u, v) = (0, 0) to (7, 5) with both left and 
right moving final steps and four right to left di- 
rection changes. 
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completely characterized by the locations of the right to left direction changes — the points 
circled in Figure 7. If there are k such direction changes, a path ending at |a; — 1) has 
amplitude 



(isin^)2'=-^(cos^)^+^-^2'=-^) = (cos^)"+"(itan^)2'=-i 
and one ending at \x) has amplitude 

(zsin^)2'=(cos^)'^+"-2'= = (cos6')"+"'(«tan6')2^ 



(11) 



(12) 



Counting the numbers of paths with k direction changes and multiplying by the amplitudes 
in (11) and (12) gives 



^t{x-l) = (cos 6*)"+^ -i tan 



k=i 



u-1 
k-l)\k 



{iiaja.9) 



2k-2 



(cos^)"+^ ■itane2Fi{l-u,-v;l\- tan^ 0) 



(13) 



and 



^t{x) = (cos6') 



S{v) + (1 - S{v)) J2 



k=l 



u 



v-1 



k \k-l 



(itan6') 



2k 



= (cos6')"+^ 5{v) - {I - 5{v))ui3n^ e 2Fi{l - u,l - v]2\ - 9) , (14) 

where the S{v) term counts the single k = path which contributes only when v = 0, and 
2F1 is the Gauss hypergeometric function. These amplitudes, together with those for an 
initially left moving path (which may be obtained from (13) and (14) by interchanging u 
and v), define K{t, x;0,0) for the quantum particle automaton. By linearity (additivity), 
therefore, they provide an exact solution for the evolution of any initial condition. 

Although Lorentz invariance is manifestly broken by the spacetime lattice, it is re- 
gained in the continuum limit: Let the lattice spacing be e, i.e., replace 9 by e9 and {u,v) 
by {u/e,v/e). Then for u,v 0, as e 0, 



and 



(j)t{x - 1) ~ ztane^o-^i(-; 1| - {uv / e^) taxi^ e9) 
r^ie9Jo(T9) 

(f)t{x) {u/e)tan^ e9oFi{-;2\ - {uv/e^)tan^ e9) 

-iue/T)MT9), 



(15) 



(16) 



where 0-^1 is obtained by taking the confluent limit of 2F1 twice [28], Ji is the ith order 
Bessel function of the flrst kind, and r := 2y/uv — — x"^ is the spacetime separation of 
(0, 0) and (t, x). The limits (15) and (16), together with their initially left moving counter- 
parts (still obtained by interchanging u and v), give exactly the continuum propagator for 
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the Dirac equation for a particle of mass 9 [29]. When ^ = 0, the contributions from (15) 
and (16) vanish and the propagator has support only on the lightcone (the 6{v) term in 
(14)). These results explain our observations about the simulations shown in Figures 1, 2 
and 3: For 9 — the propagation is at speed 1 in lattice units {i.e., is along the lightcone); 
as 9 increases the 'mass' tan^ increases and the speed decreases; when 9 = 7r/2,the 'mass' 
is infinite and there is no propagation. 

5. Two component quantum cellular automata 

In the chiral representation, the one dimensional Dirac equation describes the evolution of 
a two component complex spinor [29] . Although our interests in this paper are primarily to 
investigate possible QCA rather than to construct discrete models for fundamental physical 
processes [30], the results from the previous section indicate that in the quantum particle 
automaton it is natural to combine the amplitudes for the single particle leaving position x 
to the left and to the right into a two component field iptix) :— {(pt{x — l),(l)t{x)) for x = t 
(mod 2) and motivate consideration of two component QCA. In terms of the evolution 
rule (9) becomes 



The most general local evolution rule for a two component QCA still has the form (3), 
with (f) replaced by the two component field ip defined at all cells and the coefficients 
w now representing 2x2 matrices. The evolution rule (17) for our quantum particle 
automaton already shows that the conclusion of the No-go Lemma can be evaded in a 
two component, one dimensional QCA; the issue becomes identifying all possible local, 
homogeneous, unitary evolution matrices U as in (4). For r = 1, the most general local 
evolution rule is 




(17) 



^Jt+iix) = W-iiJtix - 1) + wo^Jtix) + w+ii/jtix + 1) 



(18) 



and the unitarity constraints are still those given by equations {5m)' 



w+^wl^ = 0, 



(19-i) 
(19o) 
(19+i) 



together with their Hermitian conjugate equations. 



Parity invariance imposes two additional constraints: 



W-i = Pw+iP 
wo = PwqP-^ 



-1 




(20) 
(21) 
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Equation (19+i) implies that at least one of W-i and tu+i is singular; by (20) both are, 
have the same eigenvalues, and hence can be simultaneously row/column reduced to, say, 

^-^=(,0 be^^ ) ^+^=(,zae^« oj- ^^2) 

That is, any parity invariant, r — 1, one dimensional, two component QCA is unitarily 
equivalent to one with w_i and in the form (22). Equation (21) implies 

Using the forms (22) and (23) in equation (19o) shows that 7 = a, (5 = and d = —cb/a. 
Finally, equation (19_i) forces 

a^ + b'' + c\l + bya^) = l 

and shows that the only nonzero solution occurs when a = (3. Thus we may factor out an 
overall phase and reparameterize to find that any nontrivial solution to equations (19^), 
(20) and (21) is unitarily equivalent to 

f ism9\ / cos ^ 

t(;_l=COSp „ „ W+i=COSp{ ■ . n n 

V COS ^ \ ismO . 

. (24) 

sin 6 —i COS 6 



wo = smp\ . 

V —I cos sm 



Our quantum particle automaton, with evolution rule given by (17), is the most general 
p = solution; the one dimensional lattice Boltzmann equation of Succi and Benzi [18] is 
unitarily equivalent to ours, as is the one dimensional version of Bialynicki-Birula's unitary 
CA for the Dirac equation [19]. As two component QCA with p = 0, each of these models 
consists of a pair of independent automata supported on the spacetime cells t + x = 
and 1 (mod 2), respectively. Setting p ^ (mod n) couples these two automata. Figure 
8 shows two simulations of a particle initially localized to be right moving from a; = 0, 
evolving with the same value of 9, but different values of p. Cell darkness is (positively) 
proportional to the probability i/j^iIj, where ip"^ is the conjugate transpose of t/j. Increasing 
the coupling p towards 7r/2 has the expected effect of decreasing the propagation speed. 

6. Quantum lattice gas automata 

Having reinterpreted our original QCA (9) as a quantum particle automaton (17), the 
lattice gas paradigm [17] suggests a different generalization than the two component QCA 
with evolution rule (18), namely to multiple particles. A LGA is a CA in which the possible 
states of the field are taken to represent occupation by ct-moving particles. The evolution 
rule consists of two stages: each particle at x jumps to x + a and then the particles now 
at X interact to possibly change their directions of motion. 



13 



From QCA to QLGA 



David A. Meyer 




From QCA to QLGA 



David A. Meyer 



In the single particle interpretation of our QCA developed in Section 4, the wave 
function 0^ = (f)t{x)\x) describes the quantum state of the system at time t, where \x) 
is the one particle state defined in (10). The scattering matrix S given in (8) may be taken 
to act on the basis states \x) rather than on the coefficient amplitudes (pt{x): 

_ / 0\x — l)t+i + i sin 6 \x)t+i t ^ x (mod 2) 
[ i sin ^ + cos ^ |a; + 1)^4.1 t = a; (mod 2). 

That is, the left (right) moving particle jumps to the left (right) and either continues in 
the same direction with amplitude cos^ or reverses direction with amplitude isin^; this is 
exactly the form of a (quantum) LGA evolution rule. 

Generalization to multiple particles has both kinematical and dynamical aspects [31]: 
The Hilbert space must be extended to have basis states | ) (denoting the config- 

uration with n particles in states xi, . . . , Xn)- The familiar restriction in classical LGA to 
occupation numbers and 1, i.e., an exclusion principle, is consistent with the fermionic 
character of the Dirac equation we found in the macroscopic limit of the quantum particle 
automaton. With this constraint the Hilbert space has dimension 2^ (if it is N for the 
corresponding QCA) since no two of the Xi in a basis state may be identical. Each particle 
still jumps at each timestep, but now there is the possibility that two particles (although 
no more than two if the exclusion principle is in effect) will jump to the same position at 
the same time. Thus the scattering matrix S must also be extended to include amplitudes 
Sij for the transitions i ^ j where 00 < i,j < 11 (in binary) and the position of the Is 
in an index indicates the occupied particle states in the pair of cells (the position) under 
consideration. Retaining parity invariance from the one particle model and imposing par- 
ticle number conservation, the most general local evolution rule is defined by the scattering 
matrix 

\ 

/ J 

which, just as in Section 2, must be unitary to ensure unitary evolution of the whole 
automaton. Dividing out an overall phase, S may be parameterized as 

/I . 

g _ ze*" sin d 

~ e*" cos d 

\ 

Thus individual particles evolve just as before except for a multiplicative phase of e*" at 
each timestep, unless two jump to the same position at the same time, in which case they 
exit multiplied by the phase e*^. This is the simplest QLGA in one dimension. 

To simulate a deterministic or probabilistic LGA, it is sufficient to store only a single 
particle configuration at each timestep and then evolve it to the next with the appropriate 



S = 



a 



b c 
c b 



le 



cos 6 
* sin^ 



(25) 
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probability (1 in the deterministic case). Each run produces a single final configuration; 
multiple runs produce the same probability distribution of final configurations as would 
computing the whole Markov process, i.e., multiplying a vector representing state proba- 
bilities by the Markov evolution matrix analogous to U. Because probabilities do not add 
in the quantum situation — there is interference in the coherent sum of amplitudes — only 
the latter procedure is viable for a QLGA. For QLGA the dimension of the Hilbert space 
is exponential in the cardinality of the lattice, so simulation is potentially exponentially 
slower than for deterministic/probabilistic LGA. With particle number conservation, how- 
ever, the problem is not completely intractable: for n fermionic particles the Hilbert (or 
Fock, as it would be called in the quantum field theory context) space has dimension 
(^) ; when n = 1 this is just the A^ amplitudes computed at each step in the one particle 
simulations of Figures 1, 2, 3 and 8. 

Figure 9 shows the next simplest situation, a QLGA with n = 2 interacting particles. 
Since one might expect that an interaction consisting only of phase multiplication (the e*^ 
term in S) would have little effect. Figure 9 compares the two particle simulation with 
particles initially at x = 4 (right moving) and x = 11 (left moving) to a one particle 
simulation with the particle initially at a; = 4 and x — 11 with equal amplitude. In each 
simulation cell darkness is (positively) proportional to the probability that a particle is 
present. The qualitative difference is immediately apparent. 

The difficulty of simulating the n particle sector of the Fock space clearly grows 
polynomially with n <^ N. But note that there is a duality between the particles and 
the 'holes' (the unoccupied states) which suggests that they should be considered to be 
antiparticles — not surprising since we saw in Section 4 that the one particle sector limits 
to the Dirac equation. In fact, the scattering matrix (25) is the Minkowski space form of 
the Boltzmann weight in the symmetric six-vertex model, which is itself a specialization 
of the Boltzmann weight for the eight-vertex model: 

d \ 



a / 

where the downward pointing arrows represent the antiparticles (the CP duals) of the 
original particles. Not only are these models exactly solvable [32], but in the appropriate 
limit (critical point) the Minkowski space symmetric six- vertex model becomes the massive 
Thirring model [33] . As usual, of course, exactly solvable does not mean that expectation 
values of all observables can be computed; even in the restricted parameter domain corre- 
sponding to probabilistic evolution only certain correlation functions have been computed 
[34] . This leaves a wide range of problems which may be most easily solved by simulation 
with this QLGA. 



y\ 

^y I a 
\.y be 
■^y e b 

■\y \ d 
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Figure 9. The QLGA with scattering matrix (25). The simulation on the left contains a single 
particle initially at a; = 4 and x = 11 with equal amplitude. The one on the right contains two 
particles, initially at a; = 4 and x = 11. In both simulations 6 = 7r/4, a = and (3 = —Stt/A. 
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7. Discussion 

Motivated by the vision of quantum computation implemented at the device level as a 
QCA, we have investigated the simplest possible models in one dimension. Although 
elementary, the No-go Lemma and its proof seem to be original. Evading its conclusion in 
order to construct a nontrivial QCA led us to the partitioning/alternating evolution rule 
which, taking the contrapositive of Henon's dictum,* we interpreted as a quantum particle 
automaton and then generalized to a QLGA. Although equivalent and similar models have 
been considered previously as regularizations of quantum field theories [18,19,33], in the 
probabilistic domain [34] , and in the context of hidden variable theories [36] , ours seems 
to be the first unitary simulation. 

Decoherence of the quantum state either in memory registers [22] or acted on by 
logical gates [23] places tight constraints on the number of computational steps which 
a quantum computer might perform reliably. QCA provide a simple model in which to 
study decoherence during time evolution. To explore the potential for quantum compu- 
tation with QCA we are currently simulating decoherence with these models, as well as 
investigating the possibility of relaxing the homogeneity condition to reflect the presence 
of local gates/devices/defects, non-periodic boundary conditions, and extension to higher 
dimensions [30]. 

These investigations are closely related to several issues in fundamental physics. From 
the assumptions of discreteness, locality, unitarity, (near) homogeneity and parity invari- 
ance we were led to a theory of chiral fermions in one dimension. In fact, the No-go Lemma 
was named to evoke the well known Nielsen-Ninomiya Theorem concerning the doubling 
of fermions on the lattice [37], although the logic here is reversed: In 1 -|- 1 dimensions 
there are several resolutions to the problem [38]; our one component partitioning QCA is 
equivalent to Casher and Susskind's [39] ; reformulating it as a two component QCA with 
two independently evolving fcrmionic automata as in Section 5 resurrects the problem; 
and generalizing to the coupled evolution rule given by (18) and (24) is analogous to Wil- 
son's solution [40]. In a more general context than regularizing continuum quantum field 
theories, the considerations involved here in imposing unitarity and locality on a discrete 
theory are also relevant to the causal set program for quantum gravity [30,41]. These con- 
nections are not surprising: the relation between physics and computation has only been 
made more intimate by the introduction of quantum mechanics. 
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"Lattice gases are at present often referred to as cellular automata. In the present note I wish to 
advance the thesis that nothing is gained by this practice, and that something is lost." [35]. 
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